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The following report summarizes the activity carried out under the 
NASA Grant entitled “ Dynamic Analysis of a Two Member 
Manipulator Arm” NAG 1-1997 from March 15, 1996 through July 6, 
1997. These activities have been completed and this document 
constitutes the final report on this project. 

This investigation was originally proposed by Dr. Elias Abu-Saba and 
he was the principal investigator when the grant was awarded. 
However, Dr. Abu-Saba retired from the North Carolina A & T State 
University, on August 15, 1995. Since Dr. Abu-Saba has retired, no 
work has been done on this project. In December, 1995, both Dr. 
Shen and McGinley expressed an interest in completing this work. 
After discussions with the technical monitor of the project. Dr. R. 
Montgomery, application procedures were started to approve 
replacement of the Principal Investigator on this NASA Grant. The 
request asked that Dr. W. Mark McGinley and Dr. Y. Shen be made 
Co-PI’s on this project. A copy of the amended proposed scope for 
this work is included in Appendix A. 

Since receiving authorization to work on this project, a computer 
model to has been developed and a number of simulations have been 
run. These simulations looked at the feasibility of using piezoelectric 
actuators to control the end affector vibrations of a two member 
manipulator arm.. Based on the results of these simulations, it 
appears that this method of damping vibrations has great potential. 



The model and results of this investigation are summarized in the 
technical paper included in Appendix B. 

During the project, a NO-COST Extension was applied for and granted 
so that the new ending date of the project was to be July 31 , 1997. 

This extension was sought to allow the completion of the attached 
paper and Dr. Shen to present these results at the Space Studies 
Institute - Conference of Space Manufacturing, May 8 - 11, 1997. 

A copy of Dr. Shen’s presentation materials are included in Appendix 
C 



APPENDIX A 


Amended Proposal 



Dynamic Analysis of a Two Member 
Manipulator Arm 
Revised Project Scope 


Dr. W. Mark McGinley and Dr. Ji Y. Shen 


The previously proposed project scope will not change.. An 
analytical model of a two member manipulator arm will be 
developed and evaluated. The results of this analysis will be 
summarized in a technical paper for presentation at a national 
conference. 

The revised budget is shown on the following page and includes 
partial summer salaries for the two investigators and travel funds. 
These travel funds are required for travel to NASA, Langley, for 
consultation with the Project Monitor and to partially cover the 
travel costs to present the technical paper at a national conference 


Proposed Project Duration 


May 5, 1996 to December 31, 1996 



Budget Dynamic Analysis of Two Member 
Manipulator Arm 

REVISED 
(JAN16, 1996) 

Item 




Principal Investigators Salaries 

Unit 

Unit cost 

Total 

Dr. McGinley Summer (per month) 

0.5 

5,497.78 

$ 2,748.89 

Dr. Shen Summer (per month) 

1 

5,222.22 

$ 5,222.22 





Sub total Salaries 



$ 7,971.11 





Fringe Benefits on Salaries (24%) 



$ 1,913.07 





Overhead (55% on salaries) 



$ 4,384.11 





T ravel 



$ 731.72 







Total 

$ 15,000.00 
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Technical Paper presented at the 
Space Studies Institute - Conference 
of Space Manufacturing, May 8 - 11, 
1997. 



END-EFFECTOR VIBRATION SUPPRESSION 
OF A FLEXIBLE MANIPULATING SYSTEM 
BY USING PIEZOELECTRIC ACTUATORS 


Ji Y. Shen, William M. McGinley and Lonnie Sharpe, Jr. 

Dept, of Architectural Engineering, College of Engineering 
North Carolina A&T State University 
Greensboro, NC 27411 


Abstract 

Attenuating start-up and 
stopping vibrations when 
maneuvering large payloads 
attached to flexible 
manipulator systems is a great 
concern for many space 
missions. To address this 
concern, it was proposed that 
the use of smart materials, and 
their applications in smart 
structures, may provide an 
effective method of control for 
aerospace structures. In this 
paper, a modified finite 
element model has been 
developed to simulate the 
performance of piezoelectric 
ceramic actuators, and was 
applied to a flexible two-arm 
manipulator system. Connected 
to a control voltage, the 
piezoelectric actuators produce 
control moments based on the 
optimal control theory. The 
computer simulation modeled 
the end-effector vibration 
suppression of the NASA 
manipulator testbed for 
berthing operations of the 
Space Shuttle to the Space 


Station. The results of the 
simulation show that the 
bonded piezoelectric actuators 
can effectively suppress follow- 
up vibrations of the end- 
effector, stimulated by some 
external disturbance. 

Introduction 

The handling of a large 
spacecraft 


Copyright of 1997 by the Space 
Studies Institute. All rights 

reserved. 

using a robotic manipulator is 
an important technology for 
future space missions. These 
operations require precision 
telerobotic maneuvering of 
large payloads using the 
Remote Manipulator System 
(RMS) of the Space Shuttle. 
During start-up and stopping, 
the direction of motion of a 
large payload is difficult to 
predict because of start-up 
transient impulses and the 
subsequent vibration in the 
system produced by the 
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flexibility of the manipulator- 
coupled system. One solution 
to this vibratory control 
problem is to conduct the 
operation slowly, in steps, and 
minimize the excitation. 
However, if objectionable 
vibrations do occur, then extra 
time is required for them to 
settle out. This solution will 
extend the time required for 
operations, and since the cost of 
orbit time is extreme high, the 
cost effectiveness of the 
mission will be reduced. In an 
effort to find more cost 
effective solutions, NASA 
continues to develop telerobotic 
technology that addresses these 
problemsj, 2)- As part of the 
effort described above, this 
investigation proposed to 
develop a methodology for 
end-effector vibration 
suppression on a flexible 
manipulator system by using 
piezoelectric actuators. 

The rapid development of 
the smart materials and their 
applications in smart structures 
shows great potential for the 
control of aerospace structures. 
Piezoelectric actuators, as a 
specific set of these smart 
materials, appear the most 
appropriate for this application. 
In general, two basic 
approaches, passive damping 
and active vibration control, 
have been studied for the 
control of smart structures^ 4 5| . 
Structures with 


bonded/embedded sensors and 
actuators made of piezoelectric 
materials are examples of 
actively controlled structures. 
In addition, the high stiffness 
of some piezoelectric materials, 
typically a ceramic, provides an 
advantage over viscoelastic 
materials as a passive damping 
mechanism. The material 
properties of the piezoelectric 
ceramic also has the advantage 
of being relatively stable with 
temperature over their 
operating range. 

A linear elastic finite 
element model for piezo-layer 
bonded beams, based on 
general finite element 
approach, had been developed 
in the authors’ previous work. 
The previous investigation 
applied a modified finite 
element model to the optimal 
controller design for vibration 
suppression of a cantilevered 
composite beam [6] . This 
approach is extended by the 
investigation described in this 
paper, and includes the 
development of a joint element 
for the end-effector vibration 
suppression of a flexible two- 
arm manipulating system. 

The system studied in 
this investigation is the NASA 
manipulator testbed for the 
research of berthing operations 
of the space shuttle to the 
space station, which consists of 
two flexible links and three 
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revolute joints. This testbed 
was assumed to be constrained 
in the horizontal plane for the 
modeling and analysis. Figure 
1 shows a schematic 
representation of the testbed. 
Each of the two links shown can 
be modeled as individual frame 
elements, some with bonded 
piezoelectric actuators. In 
response to a control voltage, 
the piezoelectric actuators can 
produce control moments and 
be used to dampen vibrations 
of the manipulator arm. 

The following sections 
describe the model 
development and summarize 
the computer simulation 
results. This simulation was 
conducted on the modeled 
manipulator system, and 
assumes that an external 
disturbance causes an initial 
deflection in one of the arms. 
This initial deflection 
stimulates a vibration of the 
manipulator system and the 
bonded piezoelectric actuators 
are used to suppress the 
follow-up vibration of the end- 


effector, using control signals 
based on the optimum control 
theory. 

Modified Finite Element 
Model for the Piezo-Layer 
Bonded Beam 

The general finite 
element model for deflection 
analysis of b e a m - 1 i k e 

structures can be characterized 
by: 

I -'/|| l '}+ 1 A'l fi’} = {/' 1 (1) 

For each beam element, the 
nodal displacement vector 
consists of the axial and lateral 
displacements and slopes at the 
two nodal points, that is, 
!v}, ={//,.«;. 0.;r tl .v +1 .6», +l } r , and the 
nodal force vector consists of 
axial forces, shears and bending 
moments at the same nodal 
points, that is, 

if), ={r.a.A/,.r . The 

stiffness matrix [/:], and the 
consistent mass matrix [m], of 
the ith beam element take the 
forms of m 
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Fig. 1 Schematic Diagram of the NASA Manipulator Testbed 
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where, /, is the length of the ith 
beam element, m i =pA i l i its 
mass, /y, its flexural rigidity. 
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If a piezo-layer is bonded 
to one surface of a beam 
element, the neutral axis of the 
composite section will change 
its location (D = depth of 
neutral axis, Figure 2), and this 
location can be determined 
based upon the force balance in 
the longitudinal direction of the 
element,, 


l [6]’ 


/) = ■ 


1 


2 1 EJi., + E h h h ] 


( 2 ) 
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Figure 2 shows a for the ith composite beam can 

schematic drawing of a piezo- be computed based upon (EJ^) p 

layer bonded beam element, and (£,/,)* through the 
where, /, is the length of the ith expression 
element, h is the width, h p and 

h b are the thickness of the +E <h ! * (3) 

piezo-actuator and the beam 

respectively, M p is the bending (Note that these I values are 

moment applied to the actuator computed about the neutral 

and M b is the bending moment ax j s 0 f the composite section.) 

applied to the beam. The 
equivalent flexural rigidity E, I , 



Fig. 2 The ith Element of a Composite Beam 


If a voltage e is applied 
to the piezoelectric actuator, an 
average normal strain, 

e p = \d /h it- , is introduced in the 

layer, where d is the electric 
charge constant of the 
piezoelectric material. This 
strain results in a normal stress 
a = e p \d /h \e , which produces a 

bending moment as given by. 
At = (bv)dv . Substitution 

of the expressions for o p and D 
into the expression for M p , we 
find that 


At =c e 

P P 


Where ’ 2 | H, K T e M 

control-moment coefficient. 


In summary, it is clear 
from the above development 
that a conventional finite 
element model can be 
formulated for the piezo-layer 
bonded beam-like structure if 
following modifications are 
made: (1) the location of the 
neutral axis of the element is 
specified by D, (Eq. 2); (2) the 
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equivalent flexural rigidity Ii t I , 
of the element accounts for the 
composite material behavior, 
(Eq. 3); (3) the ith element 
force vector {/}, includes the 


piezo-layer bending moment 
M p , that is, 

{/'}. =\v..o .. a/. + M ,r ,.o ,+At I 7 

* t * / * — i b, /< I — / - I b t¥ \ l * 


Stiffness Matrix of a 
Revolute Joint 


The function of a revolute 
joint is to connect two links of a 
kinematic assemblage. The 
connected links can have 
relative rotational motion, but 
the two nodes (say, / and ./) on 
each of the connected elements, 
respectively, remain coincident 
with each other (the 
compatibility condition). For a 
planar manipulating system, 
each node has three degrees of 
freedom, that is, the 
translational motions u and v, 
and rotational motion 6. Since 
a joint consists of two nodes / 
and ./, even though they are 
coincident, a joint will still have 
six degrees of freedom. 
Assuming that the translational 
stiffnesses are represented by 
translational spring constants 
k x and k y , and the rotational 
stiffness by rotational spring 
constant k g , the compatibility 
condition and moment 
equilibrium will produce the 
following six equations at each 
joint: (in matrix form). 
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-T 


(5) 


where, the coefficient matrix, 
designated as [&]*, is the 
stiffness matrix of a Joint R, 
and r is the joint moment. Note 
that \k |„ is singular for an 
individually joint, that is, 
D et\k\ R -0 , but it will not bring 
singularity into the global 
system, since the global system 
requires superposition to be 
applied at each joint. The 
inertia of the joint was 
assumed to be negligible and 
therefore neglected. 

Optimal Controller Design 
for End-Effector Vibration 
Suppression 

Ignoring any inherent 
passive damping effects, the 
system equation can be written 
in the state-space form, that is, 

{.V! =| ,IH .V} +!«!{«} (6) 


where, the state vector is 
defined as {.V} = {.r,,* 2 } r = {y,y} T , 
the system matrix. 


Ml = 


o / 
-M ' K 0 ’ 


the control 


influence matrix. 


\B} = 


I 

0 


0 


13 



and the control (force) vector, 

According to the optimal 
control theory,#,, the linear 
quadratic perfor-mance index, 
J, can be formulated as 

J = ^(X T OX+u T Ru)lt (7) 

The optimal controller solution 
for this linear quadratic 
regulator problem is u--R '' PX , 
where matrix ( P j can be 
obtained from the steady-state 
matrix R i c c a t i equation, 
A 1 P+PA-PBR'B 1 P=-Q. 


The solution to the state- 
space equation, (Eq. 6), can be 
obtained from a recursive 
formula provided that the 
initial conditions {A'),, are 
known, 9] , that is, 


{A'} m = +|.l| 1 (t- ,v ~/)/q ^ + ^'" 

( 8 ) 

where, the matrix exponential 
e*' is defined as 


e"' = f +A( At ) + ■ 


(Aty 


'(At) 


and consequently, 


A '(e A * - /) = l(At) + - 


I ( A/ ) 


-+ l( . A/ L+- 


However, if the optimal 
control solution has been found, 


then the state-space equation 
(Eq. 6) can be simplified as 

= (9) 

where the new system matrix 
with implementation of optimal 
control is defined as 
|“i=|.t]-j/ij[«]' , [e] r [P]. The same 
method described above for the 
full state-space equation, (Eq. 
6), can be used to solve this 
simplified version (Eq. 7). If 
the stiffness and mass matrices 
incorporate the models of the 
composite elements and joints 
described in the previous 
sections, then the resulting 
solution of Eq. 9, will model the 
behavior of the two arm 
system over time. 

Computer Simulation 

The system used to test 
the model and analysis 
techniques described 

previously was the NASA 
manipulator Testbed for the 
research of the Space Shuttle to 
the Space Station berthing 
operations. This research 
testbed is designed to model 
the berthing process, but is 
constrained to motions in the 
horizontal plane,,, Figure 1 
illustrates the principal 
components of the facility. The 
Space Station Freedom (SSF) 
Mobility Base is an existing 
Marshall Space Flight Center 
(MSFC) Vehicle that has a mass 
of 2156.4 kg, and is referred to, 
herein, as Air Bearing Vehicle 1 
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(ABV1). It represents the 
Space Station in the berthing 
operation, and is considered a 
payload on the end-effector. 
This vehicle is levitated on the 
MSFC flat floor facility using 
low flow-rate bearings. The 
other vehicle, the Space Shuttle 
(SS) Mobility Base, designated 
as Air Bearing Vehicle 2 
(ABV2), is attached to the wall 
of the flat floor facility through 
the shoulder joint, and can be 
connected to the SSF Mobility 
Base with a flexible two-arm 
manipulator system. Each of 
these arms are made of a 2.74 
m long aluminum I-beam with 
a mass of 3 7.089 kg, the 
flanges of which are 0.076 m 
by 0.0032 m and the web is 0.1 
m by 0.0032 m. There are 
three revolute joints: shoulder 

joint S , elbow joint E, and wrist 
joint W. 

Since the dimension of 
the end-effector with the 
payload can not be, in general, 
comparable with dimensions of 
the two arms, the end-effector 


will be modeled as a rigid body 
with a mass point at the wrist 
joint. The elbow joint (E Joint) 
and wrist joint (W Joint) are 
supported by air bearings. If it 
is assumed that the shoulder 
joint (S Joint) and elbow joint 
are independently driven by 
individual actuators, the control 
moments z s and x e act on the 
revolute joints S and E , 
respectively. The joint 

compliances are characterized 
by three spring constants: k x in 
x-direction, k y in y-direction, 
and k e for rotation. The 
corresponding input joint 
torques are transmitted 
through the arm linkage to the 
end -effector, where the 
resultant force and moment act 
upon the Space Station 
Freedom Mobility Base ( ABV1 ). 
For the simulation, each link 
was divided into five frame 
elements. The numbering 
system for finite elements and 
nodal points are as shown in 
Figure 3. 
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Fig. 3 A Finite Element Model of A Flexible Two-Arm Manipulator 

System 

with Piezoelectric Actuators Bonded on Some of the Elements 


The initial system 

configuration for computer 
simulation was assumed as 
follows. The up-arm (elements 
6 through 10) formed a 60" 
angle with the global X-axis, 

and the forearm (elements I 
through 5) formed a 30" angle 
with the horizontal. Because 

the two arms are not in the 
same orientation, it is 
necessary to account for the 
alignment of the two arms. The 
coordinate transformation 
matrix between the local 
element coordinate system z-y 
and the global coordinate 

system Z-Y can be expressed as 


" / in o o o o 

-/;/ / o ooo 

ooiooo 

* 7 ' “ 0 0 0 / in 0 

0 0 0 -in / 0 

0 0 0 0 0 1 

where, the direction cosines are 
/ = c o s ( x , X ) , and 
m = co s(x , Y ) .respectively . The 

element stiffness matrix in 
global coordinate system is 
then 

F'l, = V V \k ],•[£! 

where |A |, is the element 
stiffness matrix in local 
coordinate system. The same 
transformation was also 
applied to the element mass 
matrix. In addition, the 
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element nodal load vector in 
global coordinate system can be 
obtained by!/}, = | A| r where 
{/'}, is the element nodal load 
vector in local coordinate 
system. 

For aluminum, the elastic 
modules was taken 
as, a; = 7.6xio R .v/wF, and density 
as p = 2840A;s;/ni 1 . The second 
moment of area for the given /- 
beam cross-section 

I h =0.1562x10 5 /»'. For 

convenience, an equivalent 
rectangular cross section with 
height h = 0.0627 m and width 
b- 0.07 6 m is used in 
simulation, which provides the 
same value of the second 
moment of area. 

It was also assumed that 
the piezoelectric ceramic 
actuators, with thickness 
h p -0.003 w, were bonded on 
one side of the frame elements 
2, 4, 6, 8 , and / 0 . The 

piezoelectric layers were 
assumed to have the following 
properties: elastic modules 

E =63xio R A7 nr, second moment 

f> ’ 

of area / =0.224 x to "in' , the 
electric charge constant 

d = 5.3x10 m/ v . 

During the computer 
simulations both weighting 
function matrices, IQI and /Sj 
were set equal to the identity 
matrix. This assumes that each 
action will have equal weight 


on control. Further study of 
the actual values of these 
matrices may improve the 
effectiveness of this control 
strategy and should be the 
subject of further study. 

It can be assumed that an 
external disturbance causes an 
initial deflection of the forearm. 
In the simulations, the forearm 
nodes 7, 9, 11 were assumed to 
have an initial 0 .01 m 
deformation parallel to the X- 
axis and 0.1 m initial 
deformation parallel to the Y - 
axis. Nodes 8 and 10 were 
assumed to have the same 
deformations, but in opposite 
directions. These initial 
deformations would then result 
in a vibration of the 
manipulator system, when 
released. The computer 

simulation predicted the 
response of the system after 
released and under the action 
of optimal control moments 
provided by the piezoelectric 
actuators. The computational 
results show that bonded 
piezoelectric actuators can 
effectively suppress the follow- 
up vibration of the end-effector 
as shown in Figures 4 and 5, 
where Figure 4 indicates the 
decayed time history of node 
11 in X -direction, and Figure 5 
indicates the decayed time 
history of node 11 in In- 
direction. 

Concluding Remarks 
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A modified linear elastic 
finite element model, based on 
the general finite element 
method, has been developed to 
include the effects of 
piezoelectric ceramic actuators. 
This model has been applied to 
the end-effector vibration 
suppression of a flexible two- 
arm manipulator system. The 
control methodology uses an 
optimized technique that 
appears to effectively attenuate 
system vibrations and thus 
may be an effective means to 
control start-up and stopping 
vibrations present when 
maneuvering large payloads on 
space missions. 
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Computer simulation. 









The Mass Matrix and Stiffness Matrix: 




Piezoelectric Actuator’s Control Moment: 
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Optimal Controller Design 



and consequently, 



Linear Quadratic Performance Index: 




NASA Manipulator Testbed 
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Fig. 1 Schematic Diagram of the NASA Manipulator Testbed 



Computer Simulation 



The direction cosines are l=cos(x,X) and ni=cos(x,Y), respectively. 



Mechanical Properties: 




Computational Results 
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Fig. 5 The Decayed Response Time History of 
the End-Effector in the Y-Direction 



